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3 ABSTRACT 
< 

I resolved and detected an X-ray nebulosity around the recurrent nova T Pyx 
using a 98.8 ksec observation with the ACIS-S detector on-board the Chandra 

pq Observatory during the quiescent phase of the nova before its outburst in 2011. 

jn~] The nebula shows an elliptical shape with an inner semi-major axis ~ 0.45 arc 

sec and an outer semi- major axis ~ 0.9 arc sec which indicates a torus-like or 
a ring-like shell structure around the nova. There is also a (conical) elongation 
towards the southern direction of about 1.85 arc sec. This structure may be part 
of a bipolar outflow from the source/nova. The count rate of the nebulosity is 
0.0025±0.0010 c s _1 and that of the the central binary is ~ 0.003 c s _1 over the 
| 0.2-9.0 keV energy range. The best fitted spectrum of the X-ray nebula is a two- 

component plasma model (e.g., a double MEKAL) with ~ 0.6 keV and ~ 2.2 keV 

^ along with two different neutral hydrogen column densities of (0.2-0.9) x 10 22 cm -2 

and (3.0-26.0) xlO 22 cm" 2 for the two temperatures, respectively. I calculate an 
absorbed X-ray flux of (0.6-10.0) x 10~ 14 erg cm -2 s" 1 with a luminosity of (0.08- 
2.0) xlO 32 erg s _1 (at 3.5 kpc) for the X-ray nebula. The estimated shocked mass 
is <1.8x 10~ 5 M . The central source spectrum can be fitted by a single MEKAL 
model with a temperature 9.2f 54 keV yielding a luminosity of about 5.2 xlO 31 

k^j erg s -1 . The orbital period of the system is detected in the Chandra light curve. 

la 

Subject headings: X-rays: stars — radiation mechanisms: thermal — supernova 
remnants — shock waves — binaries: close — novae, cataclysmic variables — 
stars: Individual (T Pyxidis) 



1. INTRODUCTION 



Classical (CN) and recurrent (RN) nova outbursts occur as a result of thermonuclear 
runaways (i.e., explosive ignition of accreted material) on the surface of the white dwarf (WD) 
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primaries in cataclysmic variable (CV) systems ejecting material in a range 1CT 7 to 1CT 3 M 
with velocities from several hundred to several thousand kilometers per second (Shara 1989; 
Livio 1994; Starrfield 2001; Bode & Evans 2008). RN outbursts occur with intervals of several 
decades (Bode & Evans 2008). Nova outbursts show two main components of X-ray emission; 
a soft component dominating below 1 keV originating from the hot post-outburst WD and 
a hard component emitting above ~ 1 keV as a result of accretion, wind-wind and/or blast 
wave interaction (Krautter 2008). The hard X-rays are mainly caused by the shocked plasma 
emission having plasma temperatures in a range 0.1-10 keV with luminosities < afewxlO 36 
erg s _1 in the outburst stage (Balman, Krautter, Ogelman 1998, Mukai & Ishida 2001, Orio 
et al. 2001; Bode et al. 2006; Sokoloski et al. 2006; Hernanz & Sala 2002,2007; Ness et al. 
2009; Drake et al. 2009; Page et al. 2010; Vaytet et al. 2011; Orlando & Drake 2012; Nelson 
et al. 2012). There has only been one resolved and detected old CN remnant (GK Persei; 
Nova Per 1901) in the X-rays using ~ 100 ksec Chandra observation (Balman 2005; Balman 
& Ogelman 1999). In addition, some extension in the radial profiles of the X-ray emission 
was recovered using the Chandra data of the recurrent nova RS Oph, one and a half years 
after the outburst possibly associated with the infrared and radio emitting regions (Luna et 
al. 2009). Recently, Balman (2010) detected extended X-ray emission in the radial profiles 
of the recurrent nova T Pyx using the XMM-Newton EPIC pn data. 

T Pyx had five outbursts in 1890, 1902, 1920, 1944, and 1966 with an inter outburst 
time of 19±5.3 yrs (Webbink et al. 1987). A new quite delayed outburst occurred on April 
14, 2011 (Waagan et al. 2011) and was observed over the entire electro-magnetic spectrum 
including the X-rays (Kuulkers et al. 2011a, 2011b). Ground-based optical imaging of the 
shell of T Pyx shows expansion velocities of about 350-500 km s _1 (Shara et al. 1989; O'Brien 
& Cohen 1998). Hubble Space Telescope (HST; 1994-2007) observations of the shell show 
thousands of knots in Ha and [Nil] with expansion velocities of 500-715 km s~ x that have 
not decelerated and the main emission is within a radius of 5" (Shara et al. 1997; Schaefer, 
Pagnotta & Shara 2010). The joint UV, optical and the IR continuum spectrum of the 
binary system can be modeled using a blackbody of T ~ 34000 K with M ~ 10~ 8 M yr -1 , 
dominated by the accretion disk (Gilmozzi & Selvelli 2007; Selvelli et al. 2008). 



2. Observations and data reduction 

T Pyx was observed using the Chandra (Weisskopf, O'dell, & van Speybroeck 1996) 
Advanced CCD Imaging Spectrometer (ACIS; Garmire et al. 2000) for a total of ~ 98.8 ksec 
on three different pointings: 2011 January 31 (UT 16:22:36), 2011 February 2 (UT 03:05:41), 
and 2011 February 5 (UT 22:33:52) (PI=S. Balman). I used S3 (the back-illuminated CCD) 
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with the FAINT mode, and no gratings, yielding a moderate non-dispersive energy resolution. 
The ACIS has an unprecedented angular resolution of 0".49 per pixel (half-power diameter). 
The ACIS PSF has a radius of 0".418 (i.e., resolution) at 50% encircled energy which I 
exploit to recover the X-ray nebulosity around T Pyx. 

The pipeline-processed data (aspect-corrected, bias-subtracted, graded and gain-calibrated 
event lists) are used for the analysis, and acis-process-events thread is used when necessary 
with the aid of CIAO 4.3 and a suitable CALDB 4.4.2 . For further analysis, HEASOFT 
version 6.9 is utilized. In order to double check, the archived data of T Pyx in February 
2012 is also analyzed using CIAO 4.4 and CALDB 4.4.8. 

3. Analysis and Results 

3.1. The X-ray Nebula around T Pyx 

The Chandra observations of T Pyx were obtained in three sets (see section 2). The 
observations were merged with the standard procedures using a specific RA-DEC to match 
the three OBSIDs projecting the events to the nominal point of ACIS-S. Next, shifts were 
determined for the correction of the aspect solutions, updating WCS (World Coordinate 
System) and re-projecting the events matching them to the second observation (already 
projected to the nominal point of ACIS-S). At the end an aspect corrected and merged event 
file was generated. In order to create the PSF (Point Spread Function) MARX version 4.5.0 
was used to run the CHART ray tracer using the necessary off axis angles, and the source 
spectrum and exposure time. Next, the PSF rays were projected onto the detector plane 
and finally an events file was generated for the source PSF. I have analyzed the T Pyx data 
generated in 2011 (processed with pixel adjustment set to randomize) and archived in 2012 
(processed with pixel adjustment set to EDSER) and found similar results. 

The standard Chandra resolution does not properly resolve the extended emission. The 
images calculated with 0.1 pixel resolution (sub-pixel resolution) reveal the structure better. 
The images created from the merged events file and the PSF events file with the above 
pixel resolution were smoothed by 2x2 pixels using the task csmooth (using a Gaussian 
kernel) yielding a final image resolution of 0".15. Finally, the PSF image with the same 
spatial resolution was subtracted from the image of the merged events file while matching 
coordinates. The resulting subtracted images are displayed in Figure 1. The left hand 
panel shows the unsubtracted merged event file image, the middle panel shows the PSF- 
subtracted final image. In order to access the resulting nebular emission better, I also 
created an image from the brightest part of the spectrum between 1.0-2.0 keV with the same 
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sub-pixel resolution yielding a final image resolution of 0".15. The PSF-subtracted image in 
the 1.0-2.0 keV band is displayed on the right hand panel of Figure 1. The resulting images 
show extended emission with an elliptical shape. The inner semi-major axis is ~ 0".45 and 
the outer semi-major axis is ~ 0".9 . The inner semi-minor axis is ~ 0".2 and the outer 
semi-minor axis is ~ 0".4. This indicates a torus-like or ring-like structure around the nova. 
There is also a elongation towards the south extending about 1".85 from the point source. 
I note that this structure is not co-spatial with any HRMA artifacts. Moreover, I checked 
the consistency of the X-ray nebula in the individual observations by analyzing the three 
datasets and subtracting the PSF from each image. Though the S/N of the images are only 
1/3 the total image, structures resembling Figure 1 were obtained. 

I determined an approximate nebular and point source count rate by subtracting the 
PSF from the total image (merged) of the nebula within about a radius of 5". 2 . In the 
subtraction process, I used either lxl or 2x2 pixel smoothing for both the PSF and the 
total source images with the 0.1 pixel resolution, to account for the positional uncertainties 
of the photons (or anomalies in the PSF). I used a similar size photon extraction region 
obtained elsewhere in the images to account for the number of background photons. This 
yields a net count rate of 0.0025±0.0010 c s" 1 (about 40% error) for the extended emission 
and ~ 0.003 c s" 1 for the central source over the 0.2-9.0 keV energy range. 



3.2. The Chandra Spectrum of the X-ray Nebula and the Central Source 

In order to obtain a total spectrum from the three data sets of T Pyx, first, I derived 
the source+background spectra and the background spectra (using CIAO task specextract) 
for the individual observations using a photon extraction radius of 5". 2. Next, I combined 
the spectra using the CIAO task combine- spectra. The resulting total spectrum is displayed 
in the top panel of Figure 2. Since the nebular emission is less than 2" (total extend) around 
the point source a successful deconvolution of the spectrum of the nebular emission and the 
point source can not be made using the standard data extraction techniques using CIAO. 
Instead, I used the technique described below: 

The combined spectrum of T Pyx were re-binned to have 16 energy channels across 
the 1024 channels of the Chandra ACIS-S. Next, I created images for each of the given 16 
energy channel ranges using the merged events file and the PSF events file. Next, the PSF 
images created at the given energy channel ranges and smoothed by 3x3 pixels (using a 
Gaussian kernel) were removed from the images created at the same channel ranges (with 
the same smooting) using the merged events file. Resulting images were thresholded to 
have zero counts for the negative pixel values. Using these images, the number of counts 



Fig. 1. — The X-ray images of T Pyx Nebula in the 0.2 to 7.0 keV range. The resolution 
is 0."15 per pixel. North is up and West is to the right. The left hand panel is the 0.2-7.0 
keV image without subtraction of the central source PSF. The middle panel is the PSF- 
subtracted image from the left hand panel. The image on the right hand panel is the PSF 
subtracted image in the 1.0-2.0 keV range. The axes on the figures show RA (x-axis) and 
DEC (y-axis). The brightness of the three images are separately scaled. 
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in each channel range were calculated using the same extraction radius as the combined 
spectrum. The nebular spectrum was created by replacing the counts in the re-binned 
combined spectrum (see Figure 2). The same background spectrum is used for both the 
combined and the nebular spectral analysis. 

The nebular X-ray spectrum was modeled using a two-component plasma emission 
model with two different hydrogen column densities (i.e., tbabs*MEKAL+tbabs*MEKAL 
or t6a6s*PSHOCK+t6a6s*PSHOCK) since a single plasma model or a single absorption 
component yielded xl larger than 2 (tbabs-Wilms et al. 2000; MEKAL- Mewe et al. f986; 
PSHOCK-Borkowski et al. 200f). The fits with the power-law and two-temperature non- 
equilibrium ionization plasma models (e.g., NPSHOCK-Borkowski et al. 2001) resulted in 
either a large xl with the single component model fits or similar x 2 with the two-component 
model fits, yielding non-physical spectral parameters (e.g., too high temperatures or negative 
photon indices). Note that T Pyx is not a radio source in quiescence. Table 1 shows the 
spectral parameters for two-component plasma model fits with a MEKAL model (assuming 
collisional equilibrium) and a PSHOCK model (assuming non-equilibrium ionization). The 
fits yield similar parameters. The ionization timescales derived from the PSHOCK model 
indicate existence of non-equilibrium ionization plasma, but the parameter can not be con- 
strained allowing for ionization equilibrium. The spectral results (e.g., Double MEKAL) 
indicate an absorbed X-ray flux of (0.6-10.0) x 10^ 14 erg s _1 cm -2 which translates to an 
X-ray luminosity of (0.08-2) xlO 32 erg s" 1 at the source distance of 3.5 kpc (distance from 
Selvelli et al. 2008). 

In order to determine the central source spectrum, I fitted the combined spectrum with a 
composite model (tbabs*MEKAL+tbabs*MEKAL+tbabs*MEKAL). The two MEKAL mod- 
els (with the two different absorption models) represent the nebular spectrum and the last 
is used to model the central binary spectrum typical of accreeting CVs. I fixed the nebular 
spectrum parameters at their best fit values (second A^nslightly varied) and fitted the com- 
posite model to the combined spectrum. The fitted combined spectrum is displayed in the 
bottom panel of Figure 2. As a result, I derive an N H of 0.03^°' 15 xl0 22 cm 2 , a kT of 9.2z 54 
keV, and a normalization of 1.5+QgXlO~ 5 for the spectrum of the central source emission. 
This yields an X-ray flux of (1.0-5.8) x 10~ 14 erg s^ 1 cm -2 with a luminosity of (1.4-8.2) xlO 31 
erg s -1 . The count rate error of 40% in the subtraction of the nebular component is quadrat- 
ically added to the normalization, flux and luminosity. Note that the Nh value of the central 
source and the first component of the X-ray emitting nebula are in good agreement with the 
measured value of E(B-V)=0.5 (Shore et al. 2011) which yields an interstellar N H of 3xl0 21 
cm -2 . 
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Table 1: Spectral Parameters of the Nebular Spectrum of T Pyx (0.2-9.0 keV); ranges cor- 
respond to 2a errors (A% 2 = 3.84 - single parameter) for the first and 3a errors (A% 2 = 6.63 
- single parameter) for the second component; xl values of the fits are < 1.0. 



PSHOCK § x MEKAL § 2 



A^hi (xlO 22 cm" 2 ) 

kT sl a (keV) 
r! b (xlO n s cm" 3 ) 
Ki c ' d (xl0~ 5 cm" 5 ) 


n 4+0-25 
1 0+ L1 

i - u -0.4 

~3.2 
fi+ - 3 


o s +0 - 4 

u -°-0.3 

fi+ - 3 

u -°-0.3 

N/A 
fi +0 - 3 


N H2 (xlO 22 cm" 2 ) 

kT s2 a (keV) 
r 2 b (xlO n s cm" 3 ) 
K 2 c ' d (xl0~ 5 cm" 5 ) 


11.5 +22 3 5 
2.3< 
~2 
3.5+f 2 . 


1 7 0+ 9 ' 7 

1 ' - u -14.0 
o O+2.0 
z - z -1.2 

N/A 
9.8±f 5 ° 



§ x The fitted model is (tbabs*PSUOCK+tbabs*PSUOCK) using two different N H parameters. 
§2 The fitted model is (tbabs*M'EKXL+tbabs*MEKAL) using two different 7V H parameters. 

Q l keV ~ 1.2 xl0 7 K. 

b r=n t 

c The normalization constant of the MEKAL/PSHOCK plasma emission models K=(lCT 14 /47rD 2 )xEM where 
EM (Emission Measure) =J n e nn dV (integration is over the emitting volume V) . 
d & propagated count rate error of 40% is quadratically added. 
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Fig. 2. — The top panel shows the combined/total spectrum of the source plus the neb- 
ular emission in the X-rays (in red) and the nebular X-ray spectrum around T Pyx (in 
black). The bottom panel shows the Chandra ACIS-S combined spectrum fitted with 
(tbabs*MEKAL+tbabs*MEKAL+tbabs*MEKAL) model of emission. The dotted lines show 
the contribution of the three fitted MEKAL models. The lower panel shows the residuals 
between the data and the total model in standard deviations. 
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3.3. The Temporal Variations of the Central Source 

I created background-subtracted light curves with the aid of the CIAO task dmextract to 
search for any time variability. The three light curves were, then, used to generate averaged 
power spectra (PSD). I find significant modulations at 0.155±0.005 mHz which is the binary 
period of the system (and its second harmonic) above 99.9% confidence level (see Figure 3) 
where the 3a power threshold is > 77 taking into account the red noise in the PSD. The 
binary period is 1.8295(3) hrs; 0.152 mHz (Uthas et al. 2010). Note that there is a constant 
level of emission in the folded light curve at about 0.003 c s _1 (see Figure 3 right hand panel). 
Superimposed on this constant level, there is a variation between 0.003 and 0.008 c s _1 with 
a mean, net, count rate of 0.0025 c s -1 . These rates are in accordance with the rates I derive 
for the X-ray nebula and the central source, respectively. 

4. Discussion 

The imaging results indicate the existence of a torus-like or ring-like elliptical structure 
around T Pyx which most likely suffers projection effects. If one assumes a simple circular 
region viewed along a given line of sight (inclination) angle (R sin(i)= semi-minor axis) 
taking R=0".9 and semi-minor axis as 0".4, one finds an inclination angle of % < 27°. The 
inclination of the binary system is about 10° (Uthas et al. 2010). On the other hand, 
the elliptical structure could be a complicated projection effect of a conical face-on bipolar 
ejection interacting with a large spherical shell, where the cut plane will be in the form of an 
ellipse. The images of the nebula in the total and 1-2 keV band (see Figure 1) hint at such 
effects plus an elongation towards the south of 1".8-1".9. It is possible that this is a part of 
a bipolar ejection and/or a jet-like outflow from the nova. 

The structure in the X-ray nebula is consistent with the detection of a face-on bipolar 
ejection detected in the 2011 outburst of T Pyx using the near-infrared observations (Ches- 
neau et al. 2011). The orientation of the [NIIj+Ha: (HST) images, mostly composed of 
knots, also show a north to south elongation from the source position (Shara et al. 1997) 
consistent with the Chandr anebula at a larger extend. The small size of the X-ray nebula in- 
dicates that it most likely results from an interaction of the 1966 ejecta with the pre-existing 
shells/older ejecta. The size of the torus-like structure ~ 0".9 (a possible minimum due to 
projection effects) yields an expansion velocity of V exp^^OOD^kpc km s -1 . 

The extended emission in T Pyx was detected using an archival X MM -Newton EPIC 
pn observation obtained in 2006 (Balman 2010). However, the low pixel resolution (4") and 
the large PSF size (half energy radius ~ 6 // -7"; Striider et al. 2001), hindered the spatial 
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deconvolution and only distortion in the radial profiles were recovered. The extension in 
the EPIC pn image in slight excess of 15"is consistent with the HRMA+EPIC pn PSF half 
energy width of 15". Similarly, Chandra ACIS-S PSF has a half energy radius ~ 0".418 
and the Chandra nebula (+southern elongation) is in slight excess of the HEW making the 
two detections consistent. The south/southwest elongation in both detections are also in 
accordance given the differences in the instrument characteristics. 

The spectral deconvolution shows two different plasma emission components from the 
nova remnant. The calculated emission measure EM=< n e > 2 V e //, obtained from the 
normalization of the fit, yields an average electron density n e of about 32 cm" 3 and 115 cm -3 
for the colder and hotter plasma, respectively. A volume of 6.7xl0 50 cm 3 is used (consistent 
with a spherical region of 1" radius at 3.5 kpc) and a filling factor f=l (V e //=fxV) is 
assumed. If the filling factor is as small as lxlO" 2 , then the electron density can be as high 
as 320 cm" 3 and 1150 cm" 3 . The non-equilibrium ionization model (PSHOCK) indicate 
n e ~100 cm" 3 . The electron density in the X-ray nebula of GK Per was found in a range 
0.6-11.0 cm -3 for f=l (Balman 2005) and smaller than for T Pyx. 

The shocked mass in the X-ray nebula of T Pyx can be approximated as <1.8xlO~ 5 M 
assuming a fully ionized gas, and M neb ~ n e m H V e //. By comparison, Selvelli et al. (2008) 
have calculated an ejecta mass of 1O" 4 -1O" 5 M for the 1966 outburst. Shore et al. (2011) 
measured an ejecta mass of <lxlO" 5 M (f=l) for the 2011 outburst. The plasma temper- 
ature derived from the fits can be used to approximate the shock velocities using the strong 
shock relation kT s = (3/16)fim H (v s ) 2 , (T = 1.4xl0 5 v^^-i). I derive about 500-870 km 
s" 1 for the first component consistent with the minimum speed calculated from the size of the 
X-ray nebula and the expansion speeds of the 1966 outburst that are in the range 850-2000 
km s" 1 (Catchpole 1969). I suggest that the first plasma component is the forward shock 
and the second embedded component is possibly due to the reverse shock. 

Contini & Prialnik (1997) have modeled the circumstellar interaction of the T Pyx shells 
finding that as the latest shell catches the older shell already decelerated by the circumbinary 
medium a forward shock moves into the older ejecta and a reverse shock moves into the new 
ejecta. Typical densities they derive are about ~ 200 cm -3 and ~ 3000 cm" 3 for the forward 
and reverse shock zones. The model predicts a faster and hotter reverse shock than the 
forward shock. Chandr aresults are in reasonable agreement with these predictions. However, 
the ejecta have not swept up significantly more than their own mass, which makes the shell 
a Sedov remnant. There is a neutral hydrogen column density difference of about 10 times 
between the two plasma components of T Pyx which was similar to the X-ray nebula of GK 
per. This could be an absorption in a thin cold shell (e.g., between the forward and reverse 
shocks) or ionized absorption with possibly non-solar abundances (see also Balman 2005). 
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An ionized warm absorption column density, equivalent to cold absorption (by hydrogen) 
column density of lxlO 23 cm -2 , can be a factor of 10 less ~ a few xlO 22 cm -2 (see Pekon 
& Balman 2012). This two-component plasma model may, possibly, be an effect of an X-ray 
emitting non-equilibrium plasma in the nova ejecta that needs modeling. 

Selvelli et al. (2008) predict an accretion luminosity of 10 35 erg s" 1 in the X-rays for 
T Pyx consistent with the calculated accretion rates. The accretion luminosity is about 
100-500 times more than the X-ray luminosity I detect for the central source revealing a 
different origin of X-ray emission than accretion (i.e., boundary layer) (similar to the XMM- 
Newton results). I suggest that the X-rays from the point source originate from an accretion 
disk corona (i.e., or an ADAF zone) in the system. Although, the X-ray luminosity is 
consistent with the X-ray emission from stellar winds, the temperature is too high for such 
an interpretation. The orbital modulation can be caused by scattering of the X-rays from 
structures fixed in the orbital plane. 

The author thank the Chandra Observatory for performing the observations of T Pyx. 
SB thanks S. Starrfield, J.J. Drake, M. Bode, J. Krautter, M. Hernanz, J-U. Ness for crit- 
ical reading of the manuscript and valuable discussions. SB acknowledges support from 
TUBITAK, The Scientific and Technological Research Council of Turkey, through project 
108T735. 
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Fig. 3. — The left hand panel shows the PSD obtained from averaged power spectra of the 
three observations of T Pyx. The binary period and its second harmonic are labeled. The 
right hand panel is the folded average X-ray light curve of T Pyx using the detected binary 
period. 



